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We report Raman scattering measurements on iron-pnictide superconductor Ba(Fei- x Co x )2As2 
single crystals with varying cobalt x content. Upon cooling through the tetragonal-to-orthorhombic 
transition, we observe a large splitting of the in-plane phonons modes involving Fe and As 
displacements. The splitting of the in-plane phonons at the transition is strongly reduced upon 
doping and disappears for x = 0.06 qualitatively following the trend displayed by the Fe magnetic 
moment. The origin of the splitting is discussed in terms of magnetic frustration inherent to iron- 
pnictide systems and we argue that such enhanced splitting may be linked to strong spin-phonon 
coupling. 



I. INTRODUCTION 

The discovery of superconductivity in the compound 
LaOi-^F^FeAs with a maximal T c of 26 K 1 has led to 
a great excitement in the scientific community. Since 
then, an entirely new family of high-T c superconductors 
based on iron-arsenide superconducting layers has been 
discovered. Two main groups of iron-pnictides can be dis- 
tinguished : the ROMP (R = Rare-earth, O = Oxygen, 
M = transition Metal, P = Pnictogen) or 1111 family, 
and the AM2P2 (A = earth-Alcaline) or 122 family. All 
these compounds share the same building blocks, FeAs 
planes where Fe and As are in tetrahedrical coordination. 
There are strong evidences that the FeAs planes con- 
trol the transport, magnetic and superconducting prop- 
erties of these systems 5 . The undoped compounds dis- 
play magnetic order of Spin Density Wave (SDW) type 
at low temperature. They are not Mott insulators like 
undoped cuprates but rather bad metals with high resis- 
tivities. Upon doping with electrons or holes, or in some 
compounds by applying hydrostatic pressure, the mag- 
netic order and its associated structural transition are 
suppressed and a superconducting region emerges. Upon 
doping, T c can be raised until 55 K in SmOo.sFo.iFeAs 2 
even 56 K in Gd . 8 Th .2OFeAs^. 

While it is clear that the suppression of the mag- 
netic/structural transition is a requisite to the emergence 
of the superconducting order, there is no clear picture 
yet of the interplay between lattice and magnetic de- 
grees of freedom. In particular, the relationship between 
the structural and the magnetic transitions has been dis- 
cussed extensively reflecting contrasting approaches to 
the physics of the pnictides. In the strongly correlated 
or Mott-Hubbard approaches, it has been suggested that 
the magnetic and structural transitions are intimately 
connected due to competing super-exchange interactions 
in the Fe-plane which lead to magnetic frustration 4,6 . 
This frustration is believed to explain the rather low Fe 
magnetic moment. In this local moment scenario, the or- 



thorhombic distortion removes the magnetic frustration 
of the underlying tetragonal lattice. Another class of ap- 
proaches emphasizes the itinerant character of electrons 
in these systems and rather link the transition to Fermi 
surface nesting 5 . The itinerant picture is driven by the 
relative success of Density Functional Theory (DFT) cal- 
culations to predict the correct magnetic order, in stark 
contrast to the case of cuprates. It is also supported 
by the relative agreement between the calculated band- 
structure and ARPES experiments 7 . The possibility of 
strong coupling between magnetic and structural degrees 
of freedom is especially relevant since it could give clues 
to the underlying mechanism of high-temperature super- 
conductivity in the iron-pnictides. 

Here, we report a Raman scattering study of the 
lattice dynamics as a function of electron doping in 
Ba(Fei_ x Co :E )2As2. The impact of the structural and 
magnetic transition on zone-center phonons is investi- 
gated via temperature dependent measurement of the 
Raman phonon modes across the transition. Several 
anomalies are detected in the Fe or As related modes. 
The most salient one is a very large splitting of the dou- 
bly degenerate in-plane Fe-As mode which occurs at the 
tetrahedral to orthorhombic transition. The amplitude of 
the splitting, about 9 cm -1 for the undoped compound, 
is too large to be explained solely by the weak orthorhom- 
bic distortion and could be linked to strong spin-phonon 
coupling. 



II. EXPERIMENTAL DETAILS 

We have focused our study on the double lay- 
ered compound, of the 122 family Ba(Fei_ a; Co :E )2As2 
which is electron-doped with Cobalt inserted directly 
in the Fe planes. Previous experiments have deter- 
mined its phase diagram using different experimental 
techniques such as resistivity, heat capacity, magnetic 
susceptibility 8,9 , neutron scattering and X-ray diffraction 
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FIG. 1: Raman spectrum of BaFe2As2 at 305 K in the (z'x') 
configuration at A = 514.52nra. Inset : optical image of the 
undoped single crystal BaFe2As2 - geometry of the (z'x') con- 
figuration. 



measurement :: 10 ! 11 ! 12 . This compound exibits supercon- 
ductivity between x = 0.025 and x = 0.18 doping with 
a maximal T c = 24 K for x = 0.07 doping. At low 
doping, it undergoes a magnetic transition from a non- 
magnetic state to an antiferromagnetic long-range order 
associated with a Spin- Density- Wave (SDW). In the mag- 
netically ordered phase, the Fe spins are aligned antifer- 
romagnetically along a-axis and ferromagnetically along 
6-axis (a > b) in a stripe-like pattern. Very close to 
the magnetic transition, a structural one occurs from a 
tetragonal phase at high temperature to an orthorhombic 
one at low temperature. The two transitions are simulta- 
neous for undoped Ba-122 while they split upon Co dop- 
ing, the magnetic transition occurring at slightly lower 
temperatures&ii. The small orthorhombic distortion of 
the crystal structure with cooling changes the lattice pa- 
rameters of the Fe-As planes and induces different bond 
lengths along a and 6-axis 12 . 

Single crystals of Ba(Fei_ :E Co :E )2As2 were grown from 
a self-flux method by high-temperature solid-state reac- 
tions as described elsewhere 9 . Typical crystal sizes are 
2 x 2 x 0.1 mm 3 . Samples were cleaved to perform ex- 
periments on high-quality surfaces as shown in the in- 
set of Fig. [H Results on six different single crystals 
with Co doping x = 0, x = 0.02, x 0.03, x 0.04, 
x = 0.045 and x = 0.06 are reported in this study. Their 
superconducting transition temperatures were measured 
by SQUID magnetometry and by in-plane transport 9 . 
Freshly cleaved single crystals were held in a vacuum 
of 10 -6 mbar and cooled by a closed-cycle refrigerator. 
All the spectra reported here were performed using the 
A = 514.52 nm line of an Argon-Krypton laser. Typ- 
ical laser power densities focused on the sample range 
from 20 to 200 W/cm 2 (the lower power density was 
used for the lowest temperature spectra). We estimated 
the laser heating of the sample by comparing Stokes and 
anti- Stokes Raman spectra and also via the evolution of 
phonon frequencies with incident laser power at constant 



phonon frequency (z'x') (z'z') (x'z') (x'x') phonon symmetry 
124 cm^ 1 O X X X E~ g (Fe,As) 

209 cm' 1 O X O X B lg (Fe) 

264 cm" 1 | O | X | X | X | E g (Fe,As) ~ 

TABLE I: Raman peak frequencies, polarisation rules and 
phonon symmetry. 

cold finger temperature. All the temperatures reported 
here take into account the estimated laser heating. The 
scattered light was analyzed by a triple grating spectrom- 
eter ( JY - T64000) equipped with a liquid nitrogen cooled 
CCD detector. Typical Raman phonon intensities were 
very low and long acquisition times were needed (typi- 
cally 30 minutes). Incident light was polarized along z' 
direction (i.e. with a finite projection along the c-axis 
or [001] crystallographic direction) or along the x' direc- 
tion (i.e. along the [110] crystallographic direction) with 
quasi-Brewster incidence as shown in the inset of Fig. [T] 
. Scattered light polarization was collected along x' or 
z' directions. Thus four different configurations, (z'x') - 
(z'z') - (x'z') - (x'x'), were obtained by combining the dif- 
ferent incident and scattered photon polarizations. We 
note that photons polarized along the z' direction can 
probe afr-plane polarized phonons because of the finite 
projection of their polarization along the c-axis. 



III. RESULTS 

Figure [1] shows a typical Raman spectrum of BaFe2As2 
in the z'x' configuration at ambient temperature. The 
phonon modes display typical Lorentzian lineshapes with 
no detectable asymmetry. Table [I] summarizes the 
polarisation rules of the involved phonons. The un- 
doped single crystal BaFe2As2 has a tetragonal symme- 
try (IA/mmm) at room temperature. From symmetry 
considerations, one expects four Raman active phonons 
A lg (As), B lg (Fe), E g (Fe,As), E g (Fe,As). Consider- 
ing the polarisation rules given in Table [J and follow- 
ing Litvinchuk et al. 13 , we indexed the phonons peaks as 
209 cm -1 for B^ mode, which is a pure mode involving 
displacement of Fe- atoms along the c-axis, 124 cm -1 for 
the low-frequency mode and 264 cm -1 for the high- 
frequency Eg mode, which are strongly mixed modes, in- 
volving displacements of both Fe and As atoms in the 
afr-planes. We did not observe the Arsenic mode Ai g (As) 
at room temperature. 

The evolution of the Raman spectrum of undoped 
BaFe2As2 with temperature in (z'x') configuration is 
shown in Fig. O The spectra were vertically shifted for 
clarity. While the Bi g (Fe) mode simply hardens with 
cooling, a clear splitting of the low-frequency E^ phonon 
is observed between 135 K and 145 K ± 5 K. The tem- 
perature range is consistent with the magneto-structural 
transition temperature previously reported in Ba-122 12 . 
The high-frequency E^ mode is too weak to assert its 
splitting but the lowest temperature spectrum is indica- 
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FIG. 2: Top panel : evolution with temperature of the Raman 
spectrum of BaFeiAsi in the (z'x') configuration. Inset : 
zoom on the A\ g (As) mode. Bottom panel : evolution of the 
Bi^ Fe phonon linewidth with temperature. Inset : evolution 
of the Bi^ Fe phonon frequency with temperature. 



tive of a similar albeit smaller splitting. As shown in the 
inset, the K\ g (As) mode, undetectable at room temper- 
ature, gains considerably in intensity in the same tem- 
perature range (135 K - 145 K), i.e. across the phase 
transition. The considerable impact of the phase tran- 
sition on the As phonon intensity was also observed on 
Ca-122 by Choi et al^. We note that such anomaly in 
the As phonon intensity is not found in 1111 systems—^. 

Upon closer inspection, the 'Big(Fe) mode also dis- 
plays an anomaly in its linewidth at the transition as 
shown in the lower panel of Fig. [2j Contrary to Ca-122, 
the phonon frequency does not exhibit any anomaly but 
its linewidth shows a sizable decrease across the transi- 
tion. In an electron-phonon coupling model, this effect 
can be ascribed to the opening of a pseudogap in the elec- 
tronic spectrum associated with SDW formation 14 . Such 
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FIG. 3: Top panel : evolution with temperature of the Raman 
spectrum of Ba(Fei_ x Co x )2As2 in the z'x' configuration for 
x = 0, x = 0.02, x 0.03, x 0.04, x 0.045, x 0.06. 
Bottom : evolution of the E g (Fe, As) phonon frequency with 
temperature for the same dopings (on bottom). 



a pseudogap is indeed observed below 500 cm -1 in the 
far- infrared conductivity spectrum of Ba-122 by Pfuner 
et al*i£. In the following, we will argue however that spin- 
phonon coupling might play also a role in this anomaly 
and thus provide an alternative scenario. 

The most striking feature of the data presented in Fig. 
[2] is the very large splitting of the low energy phonon 
across the phase transition. At first sight, the splitting of 
the doubly degenerate in-plane phonons is consistent 
with the orthorhombic distortion that occurs at the struc- 
tural transition. At the transition, the lattice parameters 
in the afr-plane change (a > b) and the degeneracy of the 
atomic displacements involved in the E^ phonon is lifted. 
However, the amplitude of the splitting, around 7 % of 
the mode energy, is anomalously large. In fact, such a 
splitting has been already discussed and thought to be 
to small to be resolvec&i^. Indeed a simple approach 
linking the phonon frequency uo with the bond length I 
(uo 2 ~ js) 18 would yield a splitting of less than 1 cm -1 
based on crystallographic data^ and therefore cannot 
explain a phonon splitting of around 9 cm -1 . 
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The Co doping dependence of the low energy 
phonon splitting is shown in Fig. [3l where for each doping 
we report the evolution of the E g phonon frequency with 
temperature. As the Co doping increases, the splitting 
occurs at lower temperature following the trend observed 
for the structural transition. For x = 0.06, no splitting is 
observed for temperatures as low as 20 K in agreement 
with neutron and specific heat data which did not detect 
any structural phase transition for similar doping 8,11 . We 
note however that recent NMR and resistivity measure- 
ments performed on similar crystals with x = 0.06 were 
able to detect a SDW transition slightly above 30 K— 
The extracted Fe magnetic moment from NMR data was 
estimated to be very weak however : 0.05 \ib instead of 
0.9 /jlb for the undoped crystals 11 . Our data are sum- 
marized in Fig. [H where the amplitude of the splitting 
is reported as a function of both doping and temper- 
ature. Electron doping via Co doping reduces both the 
transition temperature and the amplitude of the splitting 
measured at 20 K. It decreases from 9 cm -1 at x = to 
5 cm -1 at x = 0.045. 

The amplitude of the phonon splitting, being di- 
rectly linked to the order parameter of the tetragonal- 
orthorhombic structural transition, gives us important 
information on the nature of the transition. Following 
recent neutron scattering measurements 20 on undoped 
Ba-122 showing a continuous transition for both the mag- 
netic and the structural order parameters, we have per- 
formed an order parameter analysis, Aco(T) = Auo(T = 
K) x (1 — where Auo is the phonon splitting, T s 
the transition temperature and (3 the critical exponent. 
Within our experimental accuracy, /3 values ranging 0.1 
and 0.2 were found to reproduce satisfactorily the tem- 
perature dependence of all doping levels, suggesting a 
second order-like transition at all dopings in agreement 
with neutron data of Wilson et al^ -. The temperature of 
the structural transition T s was extracted from the order 
parameter analysis and is summarized in Fig. [5l The 
evolution of T s with Co doping is in overall agreement 
with neutron scattering data of Lester et al*ii. We have 
added on the same figure the superconducting transition 
temperatures T c of the samples measured with a SQUID 
magnetometer. We note that the structural transition is 
still present in the range of Co doping where supercon- 
ductivity emerges in agreement with neutron, transport 
and NMR data2*ii*l&. 



IV. DISCUSSION 

As stated earlier, while the E^ phonon splitting is di- 
rectly linked to the structural transition, the measured 
orthorhombic distortion for x = cannot account solely 
for the large splitting amplitude observed here. This very 
large phonon splitting leads us to suspect that magnetic 
ordering has a strong impact on lattice dynamics and 
may enhance the splitting significantly. As previously 
emphasized, it is widely believed that the magnetic and 
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FIG. 4: Evolution of the amplitude of the E^ phonon splitting 
with temperature for different doping x = 0, x = 0.02, x = 
0.03, x = 0.04, x = 0.045. The dashed lines are fits using 
the order parameter power law described in the text using 
(3 — 0, 12. Inset : evolution of the amplitude of the splitting 
at low temperature (20 K) and of the magnetic moment of Fe 
with x dopin g 11 ! 12 ! 19 . 





160 




i i i i i i i i i i i i 




140 








120 




tetragonal 




100 
















80 














CD 


60 




orthorhombic \ 












40 








20 

















H 



0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 

x doping 



FIG. 5: Phase diagram of Ba(Fei- x Co x )2As2 ' temperatures 
of the structural (T a ) and the superconducting (T c ) transi- 
tions versus x doping (Co content). 



structural degrees of freedom are strongly connected in 
iron-pnictides. For example, ab-initio calculations show 
that the lattice dynamics depend strongly on both the 
Fe-spin state and the exact magnetic ordering. In par- 
ticular, better agreement is found between the phonon 
density of state observed by neutron and X-ray scatter- 
ing and ab-initio calculations when the magnetic order- 
ing is taken into account 21 . This view is reinforced by 
the observation that the doping dependence of the low 
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FIG. 6: In plane Fe-As phonon displacement patterns in 
the orthorhombic phase with the positions of the Fe and 
As atoms projected in the ab plane 13,22 . The modes arise 
from the doubly degenerate mode in the tetragonal struc- 
ture (IA/mmm) which splits into non-degenerate B23 and 
B3 5 modes in the orthorhombic structure (Fmmm). For As 
atoms, full and empty symbols are meant to differentiate be- 
tween the As atoms located above and below the Fe plane. 
The Fe-spins are drawn in black and display the stripe-like 
magnetic pattern observed by neutron scattering in the or- 
thorhombic phase 12 . The short range exchange integrals in- 
volving next-neighbors (Ji a and Jib) and nearest next neigh- 
bors (J2) Fe spins are also depicted. Note that the superex- 
change path of Ji a and Jib involves the As atoms. 



temperature amplitude of the phonon splitting follows 
a similar trend as the one reported for the Fe magnetic 
moment as shown in the inset of Fig. 2J Neutron scatter- 
ing measurements by Wilson et al*2£i also show a direct 
correlation between the temperature dependence of the 
magnetic and the structural order parameters in undoped 
BaFe 2 As 2 . 

A simple view of the magnetic order based on local- 
ized Fe spins provides an intuitive picture of the inter- 
play between magnetic and structural degrees of free- 
dom in iron-pnictides. The dominant exchange path be- 
tween Fe spins are along the a, b axis and along the di- 
agonals via As atomsi In the tetragonal phase, all the 
exchange paths are dominantly antiferromagnetic yield- 
ing magnetic frustration. The frustration is lifted in the 
orthorhombic phase where the exchange integrals along 
the a and b axis become unequivalent with a ferro and 
antiferro spin alignment along the b and a axis respec- 
tively (see Fig. [6j). Assuming a strong spin-phonon cou- 
pling, the large difference in exchange integrals along a 
and b directions^, Ji a and Ji& (see figure [6]), axis could 
in principle explain the large splitting of the in-plane 
Eg phonon 26 . A similar scenario has been proposed in 
frustrated Heisenberg antiferromagnets on a pyrochlore 
lattice where large phonon splitting have been observed 



by infrared spectroscopy across the magneto-structural 
transition 23 . In frustrated magnetic systems, the strong 
spin-phonon coupling is due to a modulation of the ex- 
change integrals by the lattice mode displacements in- 
volving magnetic ions. In these systems, the magnetic 
degrees of freedom usually drive the structural distor- 
tion, a scenario akin to the spin-Peierls transition—*^. 
In our case the non-degenerate phonon modes, as shown 
in Fig. [6l are mixed Fe-As modes and their displacement 
patterns modulate the Fe-As-Fe bond angles and thus 
all the short range exchange integrals between Fe atoms 
which occur mainly via the As orbitals. In particular, it is 
clear that the displacement pattern of each mode distorts 
differently the Fe-As-Fe angle involved in the Ji a and J15 
exchange integrals respectively. Such spin-phonon inter- 
action is expected to provide an additional contribution 
to the splitting of the E^ phonon mode and may also 
explain the dramatic enhancement of the As phonon in- 
tensity below Ts 1 ^^. 

While qualitatively consistent with our data, the valid- 
ity of such a localized spin scenario in itinerant systems 
such as iron-pnictides is however debatable 28 . In addi- 
tion in the spin-Peierls-like picture presented here, the 
phonon splitting is expected to scale with the square of 
the magnetic moment 26 in disagreement with both our 
data and neutron scattering measurements 2 ^. An itiner- 
ant approach linking the phonon displacement patterns 
to nesting properties of the Fermi surface might be more 
suitable to explain the large phonon splitting. Recent 
DFT calculations show a strong impact of Raman ac- 
tive phonons on the electronic density of state near the 
Fermi level and in particular for the E^ mode discussed 
here^ 2 .. Clearly more work is needed to understand the 
coupling between spins and phonons degrees of freedom 
in itinerant systems such as the iron-pnictides. 



V. CONCLUSION 

In conclusion, we have reported a doping depen- 
dent Raman scattering study of the lattice dynamics 
in Ba(Fei_ :E Co :E )2As2. A large splitting of the in-plane 
Fe-As phonon across the tetragonal-orthorhombic struc- 
tural transition has been observed. The phonon splitting 
is found to weaken upon Co doping and disappears for 
x = 0.06. The temperature dependence of the splitting 
below transition temperature is consistent with a contin- 
uous or second order transition. The amplitude of the 
phonon splitting cannot be accounted by the orthorhom- 
bic distortion alone and might be a fingerprint of strong 
spin-phonon coupling in iron-pnicides. 



1 Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, J. 
Am. Chem. Soc. 130, 3296 (2008). 

2 Z.A Ren, L. Wei, J. Yang, W. Yi, X.L Shen, Z.C Li, G.C 



Che, X.L Dong, L.L Sun, F. Zhou, and Z.X Zhao, Chin. 
Phys. Lett. 25, 2215 (2008). 
3 Cao Wang, Linjun Li, Shun Chi, Zengwei Zhu, Zhi Ren, 



6 



Yuke Li, Yuetao Wang, Xiao Lin, Yongkang Luo, Shuai 
Jiang, Xiangfan Xu, Guanghan Cao, and Zhu'an Xu, Eu- 
rophys. Lett. 83, 67006 (2008). 

4 Q. Si and E. Abrahams, Phys. Rev. Lett. 101, 076401 
(2008). 

5 For a review see: I. I. Mazin and J. Schmalian, Physica C 
469, 614-527 (2009). 

6 T. Yildirim, Phys. Rev. Lett. 101, 057010 (2008). 

7 H. Ding, P. Richard, K. Nakayama, T. Sugawara, T. 
Arakane, Y. Sekiba, A. Takayama, S. Souma, T. Sato, T. 
Takahashi, Z. Wang, X. Dai, Z. Fang, G. F. Chen, J. L. 
Luo and N. L. Wang, Europhys. Lett. 83, 47001 (2008). 

8 Jiun-Haw Chu, James G. Analytis, Chris Kucharczyk, and 
Ian R. Fisher, Phys. Rev, B 79, 014506 (2009). 

9 F. Rullier-Albenque, D. Colson, A. Forget and H. Alloul, 
larXiv:0903.5243l (unpublished) . 

10 D. K. Pratt, W. Tian, A. Kreyssig, J. L. Zarestky, S. 
Nandi, N. Ni, S. L. Bu d'ko, P. C. C anfield, A. I. Goldman, 
and R. J. McQueeney, arXiv:0903.2833 (unpublished). 

11 C. Lester, Jiun-Haw Chu, J. G. Analytis, S. Capelli, A. S. 
Erickson, C. L. Condron, M. F. Toney, I. R. Fisher, and 
S.M. Hayden, Phys. Rev. B 79, 144523 (2009). 

12 Q. Huang, Y. Qiu, Wei Bao, M. A. Green, J. W. Lynn, 
Y. C. Gasparovic, T. Wu, G. Wu, and X. H. Chen, Phys. 
Rev. Lett. 101, 257003 (2008). 

13 A. P. Litvinchuk, V. G. Hadjiev, M. N. Iliev, Bing Lv, A. 
M. Guloy and C. W. Chu, Phys. Rev. B 78, 060503(R) 
(2008). 

14 K.-Y. Choi, D. Wulferding, P. Lemmens, N. Ni, S. L. 
Bud'ko, and P. C. Canfield, Phys. Rev. B 78, 212503 



(2008) . 

15 V. G. Hadjiev, M. N. Iliev, K. Sasmal, Y. -Y. Sun, and C. 
W. Chu, Phys. Rev. B 77, 220505 (R) (2008). 

16 Y. Gallais, A. Sacuto, M. Cazayous, P. Cheng, L. Fang, 
and H. H. Wen, Phys. Rev. B 78, 132509 (2008). 

17 F. Pfuner, J . G. An alytis, J.-H. Chu, I.R. Fisher and L. 
Degiorgi, arXiv:081 1.2195 (unpublished). 

18 R. M. Martin, Phys. Rev. B 1, 4005 (1970). 

19 Y. Laplace, J. BoborofT, F. Rullier-Albenque, D. Colson 
and A. Forget, arXiv (2009). 

20 S. D. Wilson, Z. Yamani, C. R. Rotundu, B. Freelon, E. 
Bourret-Courchesne and R. J. Birgeneau, Phys. Rev. B 79, 
184519 (2009). 

21 T. Yildirim, Physica C 469, 425-441 (2009). 

22 M. Zbiri, H. Schober, M. R. Johnson, S. Rols, R. Mittal, Y. 
Su, M. Rotter and D. Johrendt, Phys. Rev. B 79, 064511 

(2009) . 

23 A. B. Sushkov, O. Tchernyshyov, W. RatclifT II, S. W. 
Cheong, and H. D. Drew, Phys. Rev. Lett. 94, 137202 
(2005). 

24 Y. Yamashita and K. Ueda, Phys. Rev. Lett. 85, 4960 
(2000) 

25 O. Tchernyshyov, R. Moessner, and S. L. Sondhi, Phys. 
Rev. Lett. 88, 067203 (2002). 

26 D. J. Lockwood, and M.G. Cottam, J. Appl. Phys., 64, 
5876 (1988). 

27 D. J. Lockwood, Low. Temp. Phys. 28, 505 (2002). 

28 I. I. Mazin and M. D. Johannes, Nature Phys. 5, 141 
(2009). 



